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Abstract —This paper presents a novel, model-based compressive 
antenna design method for high sensing capacity imaging appli¬ 
cations. Given a set of design constraints, the method maximizes 
the sensing capacity of the compressive antenna by varying the 
constitutive properties of scatterers distributed along the antenna. 
Preliminary 2D design results demonstrate the new method’s 
ability to produce antenna configurations with enhanced imaging 
capabilities. 

Index Terms —compressive sensing, antenna design, coded 
apertures 


I. Introduction 

Sensing systems attempt to extract as much information as 
possible about an object or region of interest by recording a set 
of independent measurements. The number of measurements, 
and the degree of their independence, determine how much 
information a sensing system can extract. Recent papers Q, 
0 have introduced the concept of a compressive reflector 
antenna for use in millimeter wave imaging applications. The 
compressive reflector antenna operates in a manner similar to 
that of the coded apertures utilized in optical imaging appli¬ 
cations 0-0 by introducing scatterers to the surface of a 
traditional reflector antenna, the compressive antenna encodes 
a pseudo-random phase front on the scattered electric held. 
By modifying the encoded wavefront from measurement to 
measurement, for example by rotating the reflector or by elec¬ 
trically changing the constitutive properties of the scatterers, 
compressive sensing techniques 0-0 can be employed with 
improved performance over the traditional reflector antenna. 
This paper describes a numerical method for optimizing the 
constitutive parameters of the scattering elements in order 
to design compressive reflector antennas with high sensing 
capacity. 

The remainder of this paper is organized as follows. Section 
[I^ shows why the sensing, or channel capacity is a reasonable 
metric to use when assessing the sensing capabilities of an 


antenna. Section [II^ describes a general compressive reflector 
antenna design approach, which optimizes the constitutive 
parameters of the scattering elements under realistic con¬ 
straints in order to maximize the sensing capacity. Section 
describes a simplifled design approach, which optimizes over 
the constitutive parameters of the reflector elements under box 
constraints. Section [V| presents preliminary 2D antenna design 
results, which demonstrate the new method’s ability to increase 
the channel capacity. Finally, the conclusions are presented in 
Section 1^ 


II. Motivation 

One of the key features of next generation sensing and 
imaging systems will be the ability to maximize the sensing 
capacity that is, the information transfer efficiency between 
the pixels in the imaging region and the data measured by 
the system. This can occur when the mutual information of 
successive measurements is as low as possible. One way to 
achieve this goal is to dynamically control the “wavefleld 
information” that is encoded in several dimensions (i.e. time, 
space, frequency, phase, polarization and novel angular mo¬ 
mentum) in any given sensing experiment — a methodology 
known as multi-dimensional codiflcation. 

The sensing matrix A of a compressive antenna working 
in a monostatic conflguration is fully characterized by its 
radiation pattern. This radiation pattern can be derived from 
the dyadic Green’s functions G(r,r',cj) of the compressive 
antenna, which can be expressed as the solution to the vector 
wave equation: 

V X — -V X G(r, r',cj) - cj^e(r, a;)G(r, r', cj) 

uj) 

= i(5 (r - r') (1) 


Using this definition for G(r, r',a;), the electric field radiated 
by a compressive antenna is given by: 

E(r, cj) = jcj J G(r, r', cj) • I(r', co’)dr' (2) 

where I(r', u;) describe the sources used to excite the reflector. 
Eq.m can be discretized into a linear system of equations, 
Eo; = G^S^, where the constants have been absorbed into 
the source term Due to the reciprocity theorem, this 

relationship can also be used to described the electric flelds 
scattered from an object of interest. In this case, the source 
distribution can be interpreted as the set of “contrast 
sources” located within the imaging region and the held 
vector are the scattered electric flelds measured by the 
receivers. The sensing capabilities of the compressive antenna 
can thus be obtained by analyzing the Green’s function matrix 
through its singular value decomposition, G^^ = UEV^, 
where U and V are orthonormal matrices and E is a rectan¬ 
gular diagonal matrix with diagonal elements > 0; the 
are known as the singular values. When a “contrast source” 
distribution v^ excites the system, an electric held is 
generated at the receivers. If the singular values are poorly 
conditioned, i.e. (Tmax/o'mm ^ 1, then it becomes difficult to 
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distinguish different source distributions from each other using 
only the measured fields. The sensing, or channel capacity of 
the Green’s function matrix can be used to compactly assess 
the imaging capabilities of the compressive antenna. For high 
Signal to Noise Ratios (SNR), the sensing capacity can be 
expressed as |T0| : 

+ (3) 

The proposed antenna design method seeks to maximize the 
sensing capacity of the Green’s function matrix by optimizing 
the constitutive parameters /i(r,cc;) and 6(r,co’) of the scat¬ 
tering elements placed upon the reflector. Since the channel 
capacity penalizes over the logarithm of the singular values, 
this approach favors systems with small condition numbers 

III. A General Design Approach 

In the optimization problem, the transmitting antenna sys¬ 
tem is described by a set of current sources located at T 
locations. Each transmitting antenna excites the M positions 
in the imaging region with stepped-frequency waveforms at K 
frequencies. The design procedure optimizes the constitutive 
properties e(r, cc) and /i(r, cc) of scattering elements located at 
N positions along the refiector. In order to allow the scattering 
elements to be dispersive, the permittivity and permeability 
of the scatterers at the k—ih frequency will be jointly repre¬ 
sented by the variable x/.. With this convention, the matrix 
Gfc (Xfe) e can be defined as the Green’s function 

matrix for sources radiating at frequency cc/c, located at the 
T transmitter positions, and evaluated at the M positions in 
the imaging region. This matrix is a nonlinear function of the 
design variables x/^. By concatenating the Green’s function 
matrices for multiple frequencies, the multi-frequency Green’s 
function matrix G(x) G can be expressed as: 

G(x) = G(xi,X 2 , .. ■,xk) 

= [Gi (xi),G2(x2),...,Gk(xx)] (4) 


the same permittivity and conductivity. As another example, 
if metamaterial scattering elements are disallowed, then the 
inequality constraint functions force the design variables to 
produce dielectric constants > 1. 

IV. A Simplified Design Approach 

This section describes how to solve a simplified version of 
Eq. In this approach, both the scatterers and the background 
medium at the scatterer locations are assumed to be non- 
dispersive and non-conductive, so that the design variables 
xi,X 2 ,... ,xx are equal and are real-valued. Moreover, the 
constraints simply restrict the electric permittivities and mag¬ 
netic permeabilities of the scatterers to lie within specified 
ranges, [eL^CR] and The simplified optimization 

problem can therefore be expressed as: 

maximize logdet (G^(x)G(x)) (6) 

subject to XL < X < Xi^ 

Eq. § can be solved efficiently using the nonlinear con¬ 
jugate gradient method GD- This method requires expres¬ 
sions for the gradient of the cost function logdetF(x) = 
logdet (G^(x)G(x)). Assuming that F(x) that is invertible, 
the partial derivatives ^logdetF(x) and are: 

Al„gdetF(x) = tr(F-‘(x)^) (7) 

( 8 ) 

A close examination of Eq. [^reveals hat the partial derivatives 
consist of the partial derivatives By defining 

Hki x) as the discretized version of the Helmholtz operator 
for frequency k, the Green’s function matrix G/c(x) can be 
expressed as: 

Gfc(x) = #H-i(x)'4' (9) 


where the vector x is the vector of concatenated design 
variables for each frequency. Assuming that M > KT, the 
channel capacity maximization problem can be expressed as 
a non-convex “max-det” problem: 

maximize logdet (G^(x)G(x)) (5) 

subject to /ig(x) < 0, g = 1, • • • , Q 

Cp(x) =0, P = ,p 

It is easy to show that maximizing the log-determinant 
of the Grammian matrix is equivalent to maximizing the 
channel capacity. Since G^^)G(x) = V(x)I]^(x)V^(x), 
logdet (G^(x)G(x)) = log (^t)’ which differs from 

the channel capacity defined in Eq. only by constants. 

The constraint functions hq (x) and Cp (x) can be non-convex 
and depend upon the specific design constraints placed on the 
dielectric scatterers. Eor example, if the scatterers are restricted 
to non-dispersive materials, then the equality constraint func¬ 
tions force the design variables xi, X 2 ,..., x^ to produce 


where $ G H/e(x) G and ^ G . 

The matrices ^ and are subsampling matrices corresponding 
to the imaging and transmitter positions respectively. Erom this 
relationship, the partial derivatives take the following 

form: 


dGkj-x.) 

dxi 




^Hfc(x) 

dxi 




( 10 ) 


The elements of the partial derivative matrix differ 

depending upon whether xi is permittivity or permeability. If 
xi is the permittivity ej at position j, then the partial derivative 
matrix takes the form: 


= Wfcdiag(l3 0 %) (11) 

where (g) is the Kronecker product and Sij G is the 
Kronecker delta function expressed as a vector, i.e. the j — th 
element of Sij equals one and all others equal zero. If xi is 
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the permeability fij at position j, then the partial derivative 
matrix takes the form: 

= -^Lediag(l3 0 Sij)L, (12) 

Ollj 11 j 

where Lc is the discretized curl operator. Computation of these 
derivatives requires K{N ^T) calls to a forward model solver 
at each iteration in order to compute the necessary Green’s 
functions. 


V. Results 


This sections presents preliminary antenna design results, 
which were generated using the simplified algorithm and 
a 2D forward model solver based on finite differences in 
the frequency domain (FDFD) (ID The design method was 
executed for two different antenna configurations, where the 
antenna operated in refiection mode and transmission mode. 
In refiection mode, dielectric scatterers are added to the 
surface of a Perfect Electric Conductor (PEC) reflector in 
order to further perturb the fields scattered by the refiector. In 
transmission mode, dielectric scatterers are placed in between 
the transmitting antennas and the imaging region in order to 
perturb the fields that otherwise radiate in the homogeneous 
background medium. 

Figures and display the design configurations for 
the refiection and transmission mode problems respectively. 
In both modes, three line source antennas, represented by 
the white circles, were used to excite the free-space imaging 
region, colored in orange. The green pixels represent the loca¬ 
tions of the scatterers to be optimized, and the red pixels in the 
refiection mode configuration represent the PEC. The antennas 
were constrained to transmit at five frequencies linearly spaced 
between 3.1GHz and 3.5GHz, and the dielectric constant of 
the scatterers was constrained to the range [1,10]; the magnetic 
permeability was restricted to /i = /io- 

Figure [^displays the optimized permittivity distribution for 
the refiection mode problem. It is important to note that the 
design problem of Eq. is non-convex, and so it is probable 
that the solution displayed in Eigure is only a locally 
optimal solution. In practice, the optimization problem can 
be solved several times using different starting points until an 
antenna design with suitable sensing capacity is found. Eigure 
2b displays the log 2 of the singular values of the Green’s 
function matrices for the original and optimized antennas. In 
this configuration, the optimized design increases the channel 
capacity by 29 bits per pixel and decreases the condition 
number by a factor of 22, from approximately 1510 to 69. 
Eigures [3a| and [3b| display the Green’s function of the middle 
antenna radiating at 3.5GHz for the original and optimized 
antennas respectively. While the phase front of the original 
antenna is fairly uniform as a function of distance from 
the transmitter, the phase front of the optimized antenna is 
noticeably perturbed. 

Eigure displays the optimized permittivity distribution 
for the transmission mode problem, and Eigure [^displays the 
log 2 of the singular values of the Green’s function matrices 
for the original and optimized antennas. In this configuration, 
the optimized design increases the channel capacity by 40 



Eig. 1: Configuration for the compressive antenna operating in (a) 
reflector mode and (b) transmission mode. Light Blue = Transmitter 
locations. Orange - Imaging region. Green = Scatterer locations. Red 
= PEC. 



(a) 



(b) 


Eig. 2: Reflection mode conflguration: (a) optimized permittivity, (b) 
log 2 of the singular values of Green’s function matrices. 


bits per pixel and decreases the condition number by a factor 
of 60, from approximately 1400 to 23. Eigures and 
display the Green’s function of the middle antenna radiating at 
3.5GHz for the original and optimized antennas respectively. 
Without the PEC reflector or any scattering dielectrics, the 
transmitting antenna radiates isotropically, and so the phase- 
front is constant as a function of distance from the antenna. In 
comparison, the phase-front of the optimized design exhibits 
a pseudo-random pattern, which conveys an increased amount 
of information compared to the isotropic phase front. 
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VI. Conclusions 

This paper describes a novel, model-based antenna design 
method for high-sensing-capacity imaging applications. By 
optimizing the channel capacity of the dyadic Green’s function 
matrix, the new approach generates antenna configurations 
with improved sensing and imaging capabilities. Preliminary 
design results using a 2D FDFD forward model for two 
antenna configurations, operating in transmission mode and 
refiection mode, demonstrate how the novel approach gener¬ 
ates antenna configurations with improved channel capacity. 
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Fig. 5: Green’s functions of the compressive antenna operating in 
transmission mode at 3.5GHz with (a) the original design and (b) 
optimized design. 


Although this paper focused on the simplified design approach, 
the theory of the general design approach allows it to be 
applied to realistic scenarios. 
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Fig. 4: Transmission mode configuration: (a) optimized permittivity, 
(b) log 2 of the singular values of Green’s function matrices. 
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